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A non-perturbative relativistic atomic theory is presented. Spin is successfully introduced in 
Schrodinger dynamics of particle motion and the resulting energy levels of Hydrogen atom are 
shown to be exactly same as that of Dirac's theory. The separation of new radial and angular 
wave equations in spherical polar coordinates is as simple as that of the Schrodinger theory. The 
solutions of the angular part of the wave equation for spin particle are exactly same as that 
of the Schrodinger theory. The spherical harmonics for spin 1/2 particle are same as in Dirac's 
theory. The radial wave function of the Schrodinger theory gets split into two spinor components. 
The new theory accounts for the energy due to spin-orbit interaction as well as for the additional 
potential energy due to spin and spin-orbit coupling. Spin angular momentum operator is as neatly 
integrated into the equation of motion as the orbital angular momentum operator. This requires 
modification to classical Laplacian operator. Consequently the Dirac matrices and the k operator 
of Dirac's theory are dispensed with. The theory also points out that the curvature of the orbit 
draws on certain amount of kinetic and potential energies affecting the momentum of electron and 
the spin-orbit interaction energy constitutes a part of this energy. It is shown that the correction 
introduced to the radial equation due to the modification of classical Laplacian is not the same as 
that of the Darwin and Pauli theroy however, the resulting energy levels of Hydrogen atom are same 
as that of Darwin, Pauli and Dirac theories. This is because the new theory is not a perturbation 
theory. Probability current density for Hydrogen atom and for free particle is discussed. The main 
focus is on the spin 1/2 bound state single electron in Coulomb potential but the theory can also be 
extended to other spin values and potentials. The Pauli equation is obtained in the non relativistic 
limit. The new wave equation is compared with the Klein-Gordon wave equation. 

PACS numbers: 03.65.Pm, 03.65.Ca, 03.65.Ge 
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1. INTRODUCTION 

It has been more than three quarters of a cen- 
tury since the first publication of Schrodinger theory Q 
and apart from Dirac's theory [_|__[ no other satisfactory 
solution to the fine structure problem of hydrogen has 
been found. Specifically, attempts at introducing spin 
in unperturbed Schrodinger formalism were unsuccess- 
ful. Inserting spin to Schrodinger formalism as a rela- 
tivistic correction, as in case of the Pauli equation with 
two spinors is a subject of perturbation theories jpjj . Sev- 
eral other relativistic wave equations dealing with various 
aspects of spin have been put forth to address large va- 
riety of problems. For example, Klein- Gordon equation 
for spin-0 particles [28|, [2!| , wave equations for describ- 
ing relativistic dynamics of a system of two interacting 
spin 1/2 particles [30l - l35| . Breit equation for two elec- 
trons [3ll [35|-[37} also called two body Dirac equation, 
generalized Breit equation for two fermions [38j, Duffin- 
Kemmer-Petiau (DKP) theory [39j - l4o| of scalar and vec- 
tor fields for describin g in teractions of relativistic spin-0 
and spin-1 bosons [46145 4j . Some of the other works in- 
clude Poincare invariant theory of classical spinning par- 
ticles [Hj], quantum mechanical embedding of spinning 
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particles [56j and spin dependent gauge transformation 
between classical and quantum mechanics [57j . All these 
theories fall under the class of perturbation theories and 
non of these theories include spin into the dynamics of 
motion. 

In this paper I have argued that the classical Lapla- 
cian operator does not provide sufficient description for 
the total momentum of the particle and needs to be re- 
vised to include the effects of spin. This is because Lapla- 
cian is more of a mathematical operator based on geome- 
try of particle motion which can very well account for the 
velocity but if there is increase or decrease of relativistic 
mass due to physical effect such as spin, the Laplacian 
cannot account for these changes. Secondly it is pointed 
out that the separation constant relating the time depen- 
dent and the time independent hamiltonian for a bound 
particle is not the same as the free particle which is com- 
monly understood to be E. For the bound particle it is 
(E + K) where K is a function of orbit curvature. This 
additional term also affects the relativistic mass with- 
out affecting the velocity. In Dirac's theory we provide 
correction to del operator from outside by introducing 
Dirac matrices, k operator and by making a proper choice 
in selecting the the radial momentum operator [T34l8j . 
In this theory, the changes have been introduced from 
within the Laplacian. In other words the spin is included 
into the dynamics of motion. The resulting energy lev- 
els of Hydrogen atom are exactly same as that of Dirac's 
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theory. 

This section provides introduction to existing 
Schrodinger and Dirac formalisms. Sec. 2 intorduccs 
unified wave equation which can be reduced to both 
Schrodinger and Dirac type formalisms. Sec. 3 dis- 
cusses energy levels and radial wave function of Hydrogen 
atom associated with improved Schrodinger formalism. 
Sec. 4 discusses current density associated with improved 
Schrodinger wave equation. Sec. 5 shows the derivation 
of the Pauli equation in the non relativistic limit. Sec. 6 
provides concluding remarks. 



A. Schrodinger wave equation 

Schrodinger Q based his wave equation on the classical 
one dimensional wave equation given by 



(hp 
dt 



= 7 



dx 2 '' 



(1.1) 



requiring that the coefficient 7 be only a constant and 
should not involve any parameters of a particular kind 
of motion of the particle. With selection of 7 = iK/2m 
and generalizing to three dimension, the non-relativistic 
Schrodinger wave equation that would satisfy the ex- 
pected wave function ip was obtained for a free particle. 



ih 



dip 



dt 2m 

Inclusion of central potential V brings it to the form 

h 2 



(1.2) 



in—— 



dt 2m 



V 2 4> + V(r,t), 



(1.3) 



This classical form of wave equation could not account 
for the fine structure of the Hydrogen atom because it 
does not include spin of the electron. Later development 
of the relativistic Schrodinger wave equation Q also had 
the same shortcoming because it could not include Pauli 
spin matrices without destroying the invariance of the 
theory. 



B. Dirac's relativistic wave equation 

Dirac [1, took another approach for addressing the 
spin problem and derived the Hydrogen energy levels that 
were in conformance with the fine structure of Hydro- 
gen. However, certain finer effects such as Lamb shift 
IH, [HI [65| are beyond the reach of Dirac's theory and 
were later addressed by quantum electrodynamics. In 
discussing the relativistic wave equations we will not deal 
with the general case of the motion of electron in the 
presence of an external electromagnetic field but we will 
start with the free particle equation and add external 
fields as required. Dirac derived his wave equation by 



linearizing the energy-momentum invariant equation in 
space derivative. 



E 2 = El + {cpf 



(1.4) 



where E = total energy of the particle and Eq = moc 2 
is the rest energy of particle. Relativistic mass is little 
used by modern physicists. Notwithstanding the modern 
usage I have used m for relativistic mass and mo for rest 
mass throughout the article. We can write Eq. (|1.4[) as, 



E = ±{(m c 2 ) 2 + c 2 (mv) 2 }^ 2 . 



(1.5) 



As usual, one can obtain Dirac's wave equation by 
putting 



{mvf = {p 2 x + P l+pl) = p 2 . 



(1.6) 



The classical relativistic hamiltonian H for a free particle 
is the right side of Eq. (|1.5|) with the positive sign. Dirac 
modified this hamiltonian to make it linear in the space 
derivative and is given by 



H = ca ■ p + j3moC 



(1.7) 



where a and (3 represent Dirac matrices. The relativistic 
wave equation is then given by 



(1.8) 



Substitution of H from Eq. (|1.7[) and replacing E by 
ihd/dt and p by — ifiV gives Dirac's wave equation jll| . 

Q \ 
ih— +ihca-V - /3m c 2 ip = 0. (1.9) 
dt J 

Dirac's equation for a central potential V is, 

ih^-+ihca ■ V -f3m c 2 - V) ib(r,t) = 0. (1.10) 
dt J 

From here on follows Dirac's theory of Hydrogen atom 
which shall not be repeated here. Introduction of radial 
momentum, radial velocity and k operators yield familiar 
Dirac energy levels for Hydrogen atom. 

In the presence of external electromagnetic field, 
Dirac's wave equation takes the form 

d \ 
zfi— — e<p — a ■ (— ificV — eA) — finite 2 j ip — 0. 

(1.11) 

where <j> and A represent electric and magnetic potentials 
respectively. This is the basis for Dirac formalism in 
short. 



2. IMPROVED SCHRODINGER FORMALISM 



We can write de Broglie wavelength in squared form 



as. 



X 2 = V 2 T 2 , (h/pf = V 2 lv 2 



(2.1) 
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where A = h/p is the associated dc Broglic wavelength [jj 
||, p momentum of particle, h Planck's Constant, V = 
c 2 /v the phase velocity, v the particle velocity, c velocity 
of light, T the period of the wave and v the frequency. 
Eqs. (|2.1[) can also be associated with the vanishing of 
the periodic invariant discussed in [58j and given by 



= A — V 2 T 2 



{hvf = V 2 p 2 



E* = V 2 P 2 , 



(2.2) 



(2.3) 



These relations can also be directly obtained from the 
fundamental relativistic mass energy equivalence expres- 
sion 



mc 2 = (E-V). 



(2.4) 



where V on the r.h.s. can represent the central coulomb 
potential. If we square and multiply and divide r.h.s. by 

v 2 , we get 



(E-Vf = ^(mv) 2 = V 



(2.5) 



In following sections we will see that in case of Hydro- 
gen atom, the energy E appearing in Eq. (|2.5[) can be 
made to represent the Dirac energy levels with a proper 
choice of quantum mechanical operator for representing 
the momentum p 2 . We can introduce the velocities cor- 
responding to the stationary states in Eq. (|2.5|) . and it 
will make no difference to this equation. Now keeping 
the velocity constant for a particular stationary state if 
we increase or decrease the relativistic mass to m + or 
m_, it will change the momentum and energy levels to 
p + , p- and E + . E_ . What is suggested here is that the 
Laplacian operator for the angular momentum has suf- 
ficient provision for altering the velocity of the particle 
but it has no provision for altering the relativistic mass of 
the particle without altering the velocity. This deficiency 
can be removed by introducing the spin operator in the 
Laplacian which I propose can alter the relativistic mass 
of the particle without changing its velocity. Therefore 
the theory based on this concept discussed in following 
sections is perfectly invariant under Lorentz group. We 
use Lorentz transformation Eq. (12. 6|) 



(E-V) = m c 2 (l-/3 2 )-i (2.6) 
with Eq. (|2.4j) to obtain the energy momentum invariant 
(E-V) 2 = [(cp) 2 + ( mo c 2 ) 2 }. (2.7) 



which is equivalent to Eq. ()1.4|) . This allows us to write 
the wave equation in the form 

Eip = ±[(c p) 2 + Kc 2 ) 2 ] 1 ' 2 ^ + Vip, (2.8) 

where ± sign indicates positive and negative energies of 
Dirac's theory. Here E represents the total orbital energy 



of the particle consisting of kinetic energy E', rest energy 
m D c 2 and the potential energy V, 

E = E' + m c 2 + V. (2.9) 

In Eq. (|2.8p we replace E by ihd/dt and p by an uncon- 
ventional quantum mechanical operator —ihVj. Assum- 
ing spherically symmetric potential, we define operator 
V 2 in spherical polar coordinates as 



d_ 

Or 



di 



(2.10) 



where L is the orbital angular momentum operator given 
by, 

d 2 



1 d ( . „d_ 

sin 000 \ de 



i 



(2.11) 



S is the spin angular momentum operator. For electron, 
S = | her where <x are the Pauli spin matrices [ll[ • If wc 
put S = in Eq. (|2.10p . we get the classical Laplacian 
V 2 of the Schrodinger theory. Operator V 2 is associated 
with the total angular momentum of the particle and is 
capable of replacing Dirac matrices. Introduction of this 
operator in Eq. (|2.8|) yields four spinors. The rectangular 
coordinates in terms of the spherical polar coordinate 
system we are using, have been defined according to the 
scheme x = r sin 6 cos (f>, y = r sin 6 sin <p, and z = r cos 8. 

The Laplacian in spherical coordinates contains the 
squared radial momentum operator p 2 which is Hermi- 
tian. However, as pointed out by several authors (l4l.l6lT- 
|64| the radial momentum operator p r of Dirac theory is 
not Hermitian. The proposed momentum operator V 2 
also contains the squared radial momentum operator p 2 
which is Hermitian. The angular momentun operator L, 
and the spin operator S contributes nothing to the the 
squared radial momentum operator p 2 . S only change 
the relativistic mass of the particle. Therefore in this 
respect the theory has an advantage over Dirac theory. 

The wave equation Eq. (|2.8[) can be written as 

- di> ±[-h 2 c 2 V 2 + (m c 2 ) 2 ] 1 / 2 iP + V^ (2.12) 



ih- 



dt 



The wave equation Eq. (|2.12[) can be considerably sim- 
plified if potential V does not depend on time. It is then 
possible to express its general solution as a sum of prod- 
ucts of functions of r and t separately. We consider a 
particular solution of Eq. (|2.12l) that can be written as 
a product ip(r,t) = z(r)f(t). A general solution can be 
written as a sum of such separated solutions. If we sub- 
stitute the above product in Eq. (|2.12[) and divide thru 
by the product, we get 



h 2 c 2 V 2 + (m c 2 ) 2 ] 1/2 z + V. (2.13) 



If we define another function u(r) such that 

1 1/2 



1 



1 



( H 



\m c 



1 



h 

m c 



1 



-V 2 u 



1/2 



(2.14) 
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then we can write Eq. f|2 . 13[) as 



* ld f 4. 
th -fdt =± 



-V 2 u+{m c 2 f 



1/2 



V. (2.15) 



Since the left side depends only on t and right side only on 
r, both side must be equal to a same separation constant 
which in this case is (E + K) where 



K = {E-V) 



-1 



nah 2 c 2 
2{E-Vf 



1/2- 



(2.16) 



K represents correction to the particle energy level due 
to the orbit curvature k. Energy levels of two particles 
having identical parameters (including velocity) , and one 
having a straight trajectory of a free particle and another 
having a curvature of the bound orbit, cannot be the 
same. This results in different relativistic masses and 
therefore different energy levels. The curvature of the 
orbit can alter the relativistic mass without changing the 
velocity, a is a function of the particle energy which we 
define as 



2{(m c 2 ) 2 -S 2 }5 



tic 



(2.17) 



In case of circular orbits k = 1/r, and for a free particle 
k = and thus K = 0. Then equation for / can be easily 



integrated to give, 

f(t) = C exp{-i(E + K)t/K), (2.18) 

where C is an arbitrary constant and the equation for u 
becomes, 

{(E - Vf - (?7i oC 2 ) 2 }it(r) = -h 2 c 2 {V 2 + Ka/2)u(r). 

(2.19) 



If we expand both sides of Eq. (|2.14[) . we find that for 
slow moving particles, z f=a u. 



1. Separation of angular and radial wave equations 



Substitution of Eq. (f2~TU|) in Eq. (|2~T9]) gives 
Eq. (|2.21[) . The radial and the angular parts are then 
separated by substituting 



u{r,e,4>) = R(r)Y{9, 



(2.20) 



in Eq. (|2.2ip and dividing thru by RY. As is expected 
in the presence of spin, the radial and the angular wave 
functions do get split into two spinor components. These 
two components are sufficient to explain both positive 
and negative enenrgies. 



u(r) 



1 d ( d 



r 2 dr V dr 



1 d ( . 8 
sin 0— 



r 2 sin 9 89 V 89 



_l d 2 

r 2 sin 2 9 8d> 



S 

rh 



u(r). 



(2.21) 



1 d ( 7 dR\ r 2 2 . 2 TrX? , ar 11 



1 d f BY 



i d 2 Y 

sin 2 9 8<t> 2 



2v\ 2 



sVy 



(2.22) 



Since the left side of Eq. (|2.22|) depends only on r and 
the right side depends only on 9 and <fi, both sides must 
be equal to a constant that we call T. Thus Eq. (|2.22[) 
gives us a radial equation 



and an angular equation 



1 



ih 



1 8 f . 8Y 
sin 9 89 V m 89 



sVy 



- TY = 0. 



1 8 2 Y 
Bin 2 9 8<t> 2 



(2.24) 



1 d 

r 2 dr 



,dR 
dr 



2 {(m D c 2 ) 2 -(E-V) 2 }R 



h 2 c 

?Lr-Lr = o, 

lr r z 



(2.23) 



The angular equation can be further separated by sub- 
stituting 



Y(9^) = 0(9)<S>(4>) 



(2.25) 
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into it and dividing by 0$. 



1 1 



9 sin OdO V hm " dS 



1 



(2.26) 



From the theory of unitary groups and infinitesimal ro- 
tations [Til ] . it follows that the spin angular momentum 
operator S 2 is strictly a constant of the motion and can 
be replaced by the number s(s + l)h 2 where s is an in- 
teger or half an odd integer and we have already defined 
r as a constant. Therefore the quantity {VT — (S/fi,)} 2 
appearing in Eq. (|2.26|) can only be a constant which we 
shall call A. Hence we end up with two equations 



d 2 $ 



de 

. smO 

sine d9 V d6 



//<!> = 0, 



A 



sin 



6 = 0. 



(2.27) 



(2.28) 



Eqs. (|2.27[) and (|2.28p have exactly the same solutions as 
that given by Schrodinger theory where v is chosen to be 
equal to square of an integer m which takes on positive 
or negative integer values or zero. Therefore, 



$m(0) = (2tt) 2 exp(im<f>). 



(2.29) 



The condition that the solution of Eq. (|2 .28[) be finite 
at cos e — ±1 limits the values of A to 1(1 + 1) where I 
is a positive integer or zero. The physically acceptable 
solutions of Eq. (|2.28[) when m = arc called Legendre 
polynomials and for m < I are called associated Legendre 
functions. Solution for A gives us the new relationship 

{Vf - (S/H)} 2 = 1(1 + 1) = L 2 /H 2 . (2.30) 

This can be written as 

{Vrti - S} 2 = L 2 . (2.31) 

This condition can be satisfied if we substitute y/Th = 
(L + S) and this will also be consistant with Eq. (|2.24[) 
which can be written as 



complete wave function which is a solution of Eq. (|2.24[) 
when r = A = 1(1 + 1) is called spherical harmonic. If 
we put r = A = 1(1 + 1) and na/2 — in the radial 
equation Eq. (|2.23p , we get the energy eigenvalues of the 
relativistic Schrodinger theory. All these solutions of the 
Schrodinger theory remain unaltered in this theory. 



2. Separation of angular wave equation for spin 1/2 
electron 

In case of hydrogen with spin 1/2 electron, the sepa- 
ration of angular wave equation [59j is carried out in the 
same manner as in Dirac's theory. This can be done by 
redefining Eq. ([2^20)) and Eq. (|232|) as 



^(r,e,4>) = e- iEt l*R n , j (r)i)f) n .(e,4>) = 

(L + s)V£,=r^. 

Here common subscripts /, 1/2 are omitted. 



(2.33) 
(2.34) 



1,(2) 



v i-l/2,m+l/2 



l + ™ + l \ 1/2 , W 
21 + 1 ) v «.Va 



; \ 1/2 

21 + 1 ) ^+1,-1/2. 



1/2 



(2.35) 



21 + 1 %m ^/ 2 



l + m + l \ 1/2 (1) 
21 + 1 ) ^"+1.-1/2- 



(2.36) 



Coefficients in Eqs. (|[2 .35|) and (|2.36[) are the Clebsch- 
Gordon coefficients |60j. Here simultaneous eigenstates 
of L 2 , S 2 , L z and S z are given by 



^m,±l/2 = Y lm(8,4>)X±i 



(2.37) 



where Yi m (9, </>) are the spherical harmonics from 
Eq. (|2.32p for S = 0, and \± are the unit spinors. 



3. Orbit curvature affects momentum 



(l + s) 2 r = th 2 y. 



(2.32) 



Furthermore, the theory of unitary groups and infinites- 
imal rotations [ill ] would require that we identify Th 2 
with eigenvalues j(j + l)h 2 of the total angular momen- 
tum operator J 2 , where j is zero or a positive integer or 
half an odd integer. This is good enough for the radial 
Eq. ([2^23]) and it is also good enough for Eq. (f2~32]> for 
integer j = I, but this will not permit seperation of vari- 
ables in Eq. (|2 . 25[) in case of hydrogen atom with spin 
1/2 electron. Therefore Eqs. (|2~2"4| thru (|2~32j) are valid 
only for integer values of j and in a special case of S = 
we get r = A = 1(1 + 1) The angular part Yi m (e, 4>) of the 



The radial wave equation Eq. (|2.23p may be rewrit- 
ten in a form that resembles the classical Schrodinger 
wave equation (|1.3|) . If we put R(r) = x( r )l r an d 
E* = E' + V with reference to Eq. (|2.9p , then the equa- 
tion for the modified radial wave function x may be writ- 
ten as 



-fr d 2 X 

— -+ 

2mo dr 2 

l(l + l)h 2 

2m,Qr 2 



V(r)- 

s(s + i)h 2 



he ( -E* 



r \ 2niQC 2 



2m$r 2 



ha.L 

2m,Qr 2 



(E* - Vf 
2mQC 2 

X, = E*x- 



(2.38) 
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Using Eq. (|2.12[) we can write 
(E - Vf = {mc 2 f = -h 2 c 2 (V 2 + na/2) + (m c 2 ) 2 , 



m 
m 



1 



-,1/2 



— (V 

m c 



kcx/2) 



(2.40) 



(2.41) 



Eq. (|2.41[) shows that the particles can have positive or 
negative relativistic mass but the rest mass is always pos- 
itive. This is the difference between particle and anti- 
particle. Since the second term in the bracket is small 
we can have 



m 
m 



n 

TTIqC 



kcx/2) 



Thus the radial motion is similar to the classical motion 
of a particle in a potential defined by the quantity within 
the large brackets. The first term is equivalent to the 
classical kinetic energy p 2 /2mo. The first, the second and 
the fifth terms are same as the classical Schrodinger the- 
ory. The fourth term is approximately equal to p 4 /8toqC 2 
which has the form of the classical relativistic mass cor- 
rection and is same as Dirac's theory. The seventh term 
includes spin-orbit coupling. The sixth term is due to the 
spin itself. The third term is associated with the curva- 
ture of the circular orbit n = 1/r which originates from 
the separation constant (E + K) given by Eq. (|2.16p . The 
magnitude of this term is dependent on the kinetic and 
potential energies in addition to the curvature. For very 
large r this term becomes negligible. Therefore there is 
a certain amount of correction to the potential energy 
when a particle travels along an orbit having large cur- 
vature. There is a corresponding change in momentum. 
This curvature term contributes only to the relativistic 
mass of the particle. The same is the case with the term 
due to spin. The classical Laplacian part of the opera- 
tor Eq. (|2.10[) contributes to both the relativistic mass 
as well as the orbital velocity of the particle. This curva- 
ture term can be convinicntly analyzed by substituting 
— E* = niQC 2 — E. The effect of introducing the separa- 
tion constant (E + K) is equivalent to adding the term We have been using two different varieties of momentum 
na/2 to the operator V 2 , provided /(*) is defined by P One associated with Vj and another with V. The 

later associates only with the relativistic mass due to 
Lorentz transformation where as the former associates 
with the relativistic mass which includes effect due to 
the Lorentz transformation plus the effect due to spin. 
For a free particle we can have na/2 = and in the 
absence of spin S = 0. This reduces Eq. (|2.41[) to 



Using Eq. (|2.10[) we can write 

9 9 1 



h 2 i 



(S 2 + 2L- S). 



Substitution of Eq. (g35) in Eq. flglgl gives 
S 2 + 2L • S - {arh 2 /2) 



m 

— = ± 

m 



1 



2 2 

p r 



2m 2 r 2 c 2 



(2.42) 



(2.43) 



(2.44) 



Eq. (I2~18l) . 

In Dirac's analysis, the third, fifth and sixth terms do 
not appear at all and the seventh spin orbit coupling term 
appear in the form of spin-orbit interaction energy asso- 
ciated with Larmor precession and Thomas precession 
• It also includes one more term as a correction to the 
potential energy and is declared difficult to demonstrate 
experimentally. 



4-. Spin-orbit interaction energy 



± 



m 



1- V - 
c 2 



-1/2 



(2.45) 



which is the familiar Lorentz transformation equation of 
special relativity. Next we substitute Eq. (|2.44|) with 
positive sign for positive energy in Eq. (|2.39[) to obtain 



In our case, the last term in brackets in Eq. (|2.38j) 
has nothing to do with the spin-orbit interaction energy 
but the last three terms together represent additional 
potential energy due to total angular momentum. Spin- 
orbit interaction energy term is hidden within the third 
term which is due to curvature effect and can be demon- 
strated as follows. It is to be noted that the parameter 
1/r in this term does not have its origin in the classical 
Laplacian. 7 in following derivation is the fine structure 
constant. 



he ( -E* 



r \2rriQC 2 



1 1 dV 



2nif\c 2 r dr 



L S 



2 2 



(arh 2 /2) 



m 2 c 2 r 2 



V 2 



2mQC^ 



he 
T 



-E* 



2rriQC i 



V 



-(E - m c 2 
moc 2 



V ( ~(mc 2 + V - m c 2 ) 



moC 



(2.39) 



V 



V27 V m o 



V 



(2.46) 

The middle term on the right has the form of the spin- 
orbit interaction energy associated with Larmor preces- 
sion and Thomas precession 0, El . Now we have suffi- 
cient information to conclude that the term na/2 appear- 
ing in Eq. (|2.16[) simulates the bound state of the elec- 
tron and provides coupling between the momentum, orbit 
curvature, spin-orbit interaction and the central coulomb 
potential. Therefore in order to deal with the free par- 
ticle, we have to put curvature k = 0, but this does not 
eliminate the spin and so it makes its presence felt in the 
presence of an external electromagnetic field as in case of 
the Pauli equation. 
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THE HYDROGEN ATOM 



Now we find a solution for F in the form 



A. Energy levels 



The second term in radial equation (|2.23|) is, 
1 



h 2 c 2 



{{mac') 2 -{E-V) 2 }R = 

(m c 2 + E-V) (m c 2 -E + V) 



(3-1) 



he 

We will substitute 

2(m oC 2 + E) 



he 



R. 



ai 



he 



and a<i 



2(m c 2 - E) 
he 



(3.2) 



Hence with respect to Eq. (|2.17[) . a 2 = ai«2- The at- 
tractive coulomb interaction between an atomic nucleus 
of charge +Ze and an electron of charge — e is repre- 
sented by the potential energy V(r) — —Ze 2 k' jr where 
k' is Coulomb's constant. In case of hydrogen atom, 
Z = 1. We introduce these hydrogen parameters into 
Eq. (|3.1| with the fine structure constant 7 defined as 
7 = (e 2 k'Z)/{hc). 



h 2 6 



■{(m c 2 ) 2 



(E- Vf}R 

7 f a 2 - a 1 
2a 2 



(3.3) 



R 



We will rewrite the radial equation (|2.23|) in dimen- 
sionless form by introducing a unitless independent vari- 
able p = ar which we have already defined earlier with 
Eq. ([2~T7| . Substitution of Eqs. (031) and f2~3"2]) gives, 



1 d 



p 2 dp \ dp 



,dR 



P 



5- 



7 l 0i2 — ot\ 
P 



2a 
j(j + 1) 



R 



(3.4) 



R = 0, 



As far as the leading terms are concerend, for suffi- 
ciently large p it is apparent that R(p) = p n e ziz 2 p satisfies 
Eq. ([3.4)1 when n has any finite value. This suggests that 
we look for an exact solution of Eq. (|3.4|) of the form 



R(p) = F(p)e-*t 



(3.5) 



where F(p) is a polynomial of finite order in p. Substi- 
tution of Eq. p.5[) into Eq. Q3.4p gives equation for F{p) 
as 



F" 



-IF' 



A - (1/2) j(j + 1) - 7 2 



where we have substituted 



A = -7 



«2 - "1 

2a 



F = 0, 
(3.6) 

(3.7) 



F{p) = p s (a + a lP + a 2 p 2 + •■■)= p s L{p), (3.8) 

a ^ 0, s > 0. (3.9) 



Substitution of Eq. ([3.8)1 into Eq. ()3.6|) gives us the equa- 
tion for L. 



p 2 L" + p{2(s + 1) - p}L' + [p{\ - s - (1/2)} 



+ s 2 + s + 7 2 



■j(j + l)]i = 0. 



(3.10) 



If we set p = in Eq. ([3~TU[) . it follows from Eqs. ([3~g)l 
and BUI) that 



S 2 + s + 7 2 -j(j + l) = 0. (3.11) 
This quadratic equation in s has two solutions, 
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(3.12) 



The boundary condition that R(p) be finite at p = 
requires that we choose upper sign for s. It is to be noted 
here that both, Schrodinger's relativistic theory as well 
as Dirac's theory allow value of s which is very slightly 
less than the permissible value. This problem does not 
exist in this theory. The smallest value in this theory is 
~ (1 - 7 2 )/2, well within the range s > 0. With this, 
Eq. p. 101) reduces to 



pL" + {2(s + 1) - p}L' + {A 



(1/2)}L = 0. (3.13) 



Equation ([3.13p can be solved by substituting Eq. 
The recursion relation between the coefficients of succes- 
sive terms of the series is observed to be 



(is + s + (1/2) - A) 
v(v + l) + 2(v+l)(s + l) 



(3.14) 



If the series does not terminate, its dominant asymptotic 
behavior when p — > 00 can be inferred from the coeffi- 
cients of its high terms: 



1 



(3.15) 



This means that the series has the asymptotic form e p 
and regular solution is obtained only if it terminates. 
Suppose that this occurs at v = n' , so that a n >+i = 0. 
Then Eq. p~TI| give the relation 



s + (l/2)-A = 0, n' = 0,l,2, 



(3.16) 



Substituting for s and A from Eqs. ([3.12)1 and ([3.7[1 re- 
spectively, we get the energy levels of hydrogen atom 



E = TTIqC 



(3.17) 
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where the radial quantum number n' is related to the 
total quantum number n by the expression 



n' + (j + (1/2)). 



(3.18) 



These are exactly the same Dirac energy levels having the 
same total spread in energy of the fine-structure levels for 
a given n. 



B. Radial wave function 



It is clear that Eq. (|3 . 1 3|) has the form of associated 
Laguerre differential equation [ill ] , 



pLP»+(p+l-p)LP' + (q-p)Lr = 0. 



(3.19) 



The associated Laguerre polynomials LP q can be con- 
structed according to formula 



Lq (P) = E(- 1 )" +P 



[ 9 !]V 



k=0 



(q - p - k)l(p + k)lk\ 



(3.20) 



In case of Eq. (|3.13|) , we have p = 2s + 1 , (q — p) = A — s — 
(1/2) and q = A+s+(l/2). Comparision with Eq. (|3T6| 
shows that (q—p) are integers but p and q are not integers. 
Therefore it is possible to solve Eq. Q3.13P using Laguerre 
polynomials only if we introduce the approximations p « 
2j + 1 and q ps A — s + 2j + (1/2) by ignoring j 2 term in 
s. Hence we get, 



(p) 



T 2j+1 

^X-s+2j+(l/2) 
X-s-(l/2) 

(-_ 1 ^+2j+i 

k=0 



[(A - s + 2j 



(1/2))!]V 



(A - s - (1/2) - k)\(2j + 1 + k)\k\ 

(3.21) 



for j + (1/2) in Eq. (|3.4[) we can write 



1 d 



p 2 dp \ dp 



,dR 



A 

P 
1 

Yp 



k(k + 1) - y 



R 



(1/4) 



(3.23) 



R = 0, 



If one neglects the last term, the remaining two terms 
constitute the ordinary relativistic Schrodinger equation 
for particle that does not have spin. Therefore the third 
term represents the correction introduced by this theory. 
If we replace (fc + 1) by — k in Eq. (|3.23[) . we do not alter 
the first two terms but the third term will change to 



1 fc+(3/4) \ 
2p p 2 ) 



R, 



(3.24) 



which is the second possible correction to the ordinary 
relativistic Schrodinger equation. Now, even by trans- 
forming the wave function R, it is not possible to relate 
these two correction terms with that of Pauli-Darwin and 
Dirac's theories which are given by 



e 2 k' 



2nioc 2 r 



and 



e 2 k' 



2moc 2 r 3 



(fc + 1) 



(3.25) 



At the same time the energy levels of the Hydrogen atom 
are exactly same as that of Dirac's theroy. The principle 
reason being that in present theory, the spin is included 
in the dynamics of motion where as in Pauli-Darwin and 
Dirac theories it is introduced as a perturbation. For the 
same reason the spin-orbit interaction energy appears in 
Dirac's theory in a little different way than in this theory. 

Next we compare wave equation of this theory with 
the Klein-Gordon wave equation which is the square of 
the Dirac equation and can be written as 



(n 2 + p 2 )^j = o 



where 



/' 



— , (3.26) 



This will yield two associated Laguerre polynomials L p q 
corresponding to two values of A associated with two 
spinors which are applicable to positive energy as well as 
the negative energy solutions. The resulting approximate 
radial wave function is of the form e~^ p p s L 2 J +1 . „ . , 
and the normalization constant may be found by using 
the generating function to evaluate the integral 



o 2s [L 2] 



A-s+2j + (l/2) 



(p)] 2 P 2 dp. 



(3.22) 



C. Comparision with previous theories 



with previous theories of Schrodinger |7|, Darwin-Pauli 



We shall now compare the radial equation Eq. p.4[) 
with previous theories of Schrodinger [7[, 
[H Gj] and Dirac @]. Into ger eigenvalues of the k op- 
erator in Dirac's theory are same as the orbital angular 
momentum k in Pauli-Darwin theory. If we substitute k 



1 d 2 
□ 2 = — — - V 2 

c? dt 2 



d'Alembert operator. (3.27) 



With respect to the wave equation Eq. (|2.12j) . for a free 
particle we have V = and k = 0. Hence 



where the first two terms denote the modified d'Alembert 
operator D 2 which includes the spin operator and hence 
can account for all spin values. 



(□ 2 + m 2 ) ^ = 0. 



(3.29) 



The Klein-Gordon wave equation (|3.26[) cannot account 
for spin 1/2 particles, it is beleived to account only for 
scalar spin particles. The wave equation of this theory 
Eq. (|3.29[) can account for all spin values: 0, 1/2, 1 and 
2. The spin is just a quantum number. 
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4. CURRENT DENSITY 

The total angular momentum operator V 2 given by 
Eq. (|2.10l) can be written as 



y2^ V 2__gl_ 2L - S . 

3 h 2 r 2 h 2 r 2 



(4.1) 



Substitution of Eq. (|4~Tj) in Eq. (|2~19|) gives 

{(E - V) 2 - {m c 2 ) 2 }u(r) = -h 2 c 2 
2L- S 



h 2 r 2 



(4.2) 



u(r). 



{E + K)z(y) =ihc 



2L • S 



h 2 r 2 h r 



(m c 2 f 
H 2 c 2 



(4.3) 



(E + K)z{r) = 

_ h 2 v 2 



z(r) + Vz(y). 



2L • S 



z(r) + Vz(r). 

(4.4) 

Expanding this to first order accuracy we can write 



it) 



djp_ 

at 



^2 r 



2mr 



2L • S S 2 



h 2 r 2 H 2 r 2 



ip + Viji, (4.5) 



where ih(d/dt) = E* = (E + K — m c 2 ). In order to 
evaluate the spin-orbit coupling term, we make use of 
the following identity applicable to Pauli matrices. 



(er • B)(er • C) = (B • C) + i<T ■ (B X C). 



(4.6) 



Eq. (|4.6[) can be verified using similar identity applicable 
to Dirac matrices. 



(a ■ B)(a ■ C) = (B ■ C) + ia' ■ (B x C). 



(4.7) 



If we replace vectors B and C in Eq. (|4.6[) by vectors p 
and r, we get 

(o--p)(<7T) = (p-r)+Mr-(pxr). (4.8) 
icr ■ (p X r) = (iW • r) - (cr • ihV){a ■ r). (4.9) 



a ■ (p X r) = H [(V • r) - (tr • V)(er • r)] 



(4.10) 



Therefore the spin-orbit coupling term in Eq. (|4.5p can 
be written as 



2L • S _ha- -(pxr) 

= ft 2 [(V-r)-(,T.V)(,T.r)] 

/i 2 r 2 



(4.11) 



Substitution of Eq. (|4.1ip in Eq. (|4.5[) gives 



<9i 2mo 



V»-(V^) + (V. ff )(^«r) 



J2 1 



h 2 r 2 



(4.12) 



V> + Vip. 



The hermitian adjoint of Eq. f)4. 12[) is written as 
ih 



2mo 
S 2 



(V .-)-(V.<r)(-.<x) 



(4.13) 



We multiply Eq. P~T2"|) on the left by ■0 t and Eq. P~T5j) 
on the right by ip and add one to another. 

|(^) = |(^ 1+ ^ 2 ) = ^ = 

^v 2 ^ - ^vV - <jV- ■ (W) + v>- ■ (vv> f ) 

2mo [ I r r ' 



(4.14) 



If we substitute cr 2 = 3 in Eq. (|4.14j) , it reduces to 

dw ih ^ 

dt 2mn 



r r 



Eq. (|4. 1 5[) has the form of continuity equation 
dw _ . 



(4.15) 



(4.16) 



provided we define the probability current density vector 
j as 



J 



2imo 



r r 



(4.17) 



From Eqs. (|4.15|) and (|4. 1T[) it is clear that in case of 
a bound electron, the spin-orbit coupling term does con- 
tribute to the probability density and the current density. 
For a free particle we can assume a very large value for 
r in these two equations and the last two terms become 
negligible, hence we get for two spinors 



dw ih ^ 

dt 2m 



(4.18) 



2imo 

h 

2imo 



(4.19) 
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This is same as the first term of the current density for 
Pauli equation. In comparision the free particle current 
density for spinlcss particle in Schrodinger theory is 



S(r,t) 



2imn 



(4.20) 



and Dirac's theory with four spinors yields 

S(r,t) = ajj^otip, (4.21) 



In case of Schrodinger current density Eq. (|4.20p . the 
classical results for Hydrogen atom indicate nil value [2(| 
when the quantum number m = 0. In case of Dirac 
current density Eq. (|4.21|) , this value is not nil and corre- 
sponds to a circular rotation of probability [27} ■ A good 
approximation of the Dirac current is obtained when 
Gordon current [21( Si = (h/2mo)rot(tp*eT^A is added 
to Eq. (|4.20[) in the Pauli approximation [22], [H[ and 
Si = (1 /too) VP x S in the Schrodinger approximation 
pil [2a |. where P is the probability density and S the 
constant spin vector (2(| . 



5. THE PAULI EQUATION 

In the non-relativistic limit, for the motion of electron 
in the presence of an external electromagnetic field, the 
present theory can be shown to reduce to the Pauli equa- 
tion In the non-relativistic limit, we have na/2 = and 
S = 0, and the equation of motion is given by Eq. (|2.19[) 
where V = ecf>. We use the following standard quantiza- 
tion procedure in the presence of an external electromag- 
netic field. V = (81,82,83), 



lh 8l e " 



-ihV - -A 

c 



(5.1) 



where e is the charge on the particle and c is the speed of 
light. The electric and magnetic field strengths arc given 
in terms of the respective potentials </> and A by 



E 



1 8 A 
c~8t 



V4>, 



H = V x A. 



(5.2) 



Therefore Eq. (j2~19|) together with Eq. (j2~T0| reduces to 



{{E 



(m D c 2 ) 2 }?/> = (cp- eA) 2 V>, 



(5.3) 



Here we make use of the identity Eq. (|4.6[) and replace 
both vectors B and C with the vector (cp — eA). This 
gives 



[a- (cp-eA)] 2 = (cp-eA) 2 

+ icr ■ [(cp - eA) x (cp — eA)]. 

Using Eq. (|5.2j) we can derive the relation 

(cp — eA) x (cp — eA) = — ce(A x p + p x A) 
= iehcV x A = iehcH. 



(5.4) 



(5.5) 



Substitution of Eq. (|5.5|) in Eq. (|5.4[) gives 

(cp - eA) 2 = [er • (cp - eA)] 2 + ehc a ■ H. (5.6) 

In case of the motion of electron in an external electro- 
magnetic field, Eq. (|2.9I) gets modified to 



E = E' 



m c 



ecj). 



(5.7) 



If we put E* = E' + ecj> and assume that E* and ecj) are 
small in comparision with m c 2 , we can then make the 
replacement 

[E - ecj)) 2 - (m c 2 ) 2 ps 2m c 2 (E* - ecj)), (5.8) 
in Eq. (|5.3p and substitute Eq. (|5.6p to obtain 



E*if> 



1 



2mo 



^A 

c 



ecf> - 



eh 
2toqC 



er H 



where tp 



01 
02 



(5.9) 



(5.10) 



Eq. (|5.9[) is a system of two coupled differential equations 
for ipi and "02, describing electrons with z component of 
their spin up or down respectively. Therefore we can 
replace er 2 by a 2 z = 1 and in the non-relativistic limit, 
Eq. (|5.9p reduces to the Pauli equation 



L (p-" a ) : 

Tin V C / 



2toq 



ecj) ■ 



2moc 



H 



0, 

(5.11) 



where the H term represents the magnetic moment of an 
electron. 



6. CONCLUSION 

Spin is successfully introduced in Schrodinger formal- 
ism. The resulting energy levels of Hydrogen atom in 
this theory are exactly same as that of Dirac's theory. 
The separation of the new wave equation in spherical 
polar coordinates is as simple as that of the Schrodinger 
theory and the solution of the angular part of the wave 
equation for spin 1/2 electron is same as in Dirac's 
theory. The separation constant relating the time 
dependent and the time independent hamiltonian for a 
bound particle is not the same as the free particle but 
has an additional term which accounts for the energy 
contribution in the formation of the orbit curvature. 
The radial wave function of the Schrodinger theory 
gets split into two spinor components. These two 
components are sufficient to explain both positive and 
negative enenrgies. The new theory accounts for the 
energy due to spin-orbit interaction as well as for the 
additional potential energy due to spin and spin-orbit 
coupling. Spin angular momentum operator S is as 
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neatly integrated into the equation of motion as the 
orbital angular momentum operator L. Consequently 
the Dirac matrices and the k operator of Dirac's theory 
are dispensed with. The theory also points out that 
the curvature of the orbit draws on certain amount of 
kinetic and potential energies affecting the momentum of 
electron and the spin-orbit interaction energy constitutes 
a part of this energy. A new relationship between the 
rclativistic mass and the rest mass is discussed which 
reduces to the Lorentz transformation equation in the 
absence of the spin. It is shown that the correction 
introduced to the radial equation due to the modification 
of classical Laplacian is not the same as that of the 
Darwin and Pauli theroy however, the resulting energy 
levels of Hydrogen atom are same as that of Darwin, 
Pauli and Dirac theories. This is because the new theory 
is not a perturbation theory. Both the spin and the 
curvature effect are assumed to affect only the relativistic 
mass of the particle and the Laplacian operator part 
of Eq. (|2.10[) contributes to both the relativistic mass 
as well as the orbital velocity of the particle. Dirac's 
theory cannot detect the curvature effect because the 
spin in that theory is not the part of the dynamics of 



motion but it gets introduced as a perturbation just 
like Darwin and Pauli theory. Also the selection of the 
radial momentum operator is somewhat arbitrary in 
that theory [3411 and not hermitian as pointed out by 
several authors. The probability density and the current 
density for the hydrogen atom have been calculated. 
The probability density and the current density for a 
free particle is comparable to that of the Pauli equation 
for the motion of electron in an electromagnetic field 
with spin included. In the presence of the external 
electromagnetic field, the theory reduces to the Pauli 
equation in the non relativistic limit. Comparision of 
the new wave equation with the Klein-Gordon equation 
shows why the later cannot account for spin values other 
than zero. 
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